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1.  Introduction 


The  electrical  characteristics  of  an  arc  that  is  produced  by  a  high-current  pulse  have 
been  studied  for  spark  gap  switches,  lightning  and  thunder,  exploding  thin 
conductors,  and  ignition  of  a  fuel/air  mixture  in  an  internal  combustion  engine.  All 
of  these  arcs  have  a  circular  cylindrical  channel  of  an  expanding  gas  that  is  nearly 
completely  ionized  at  a  high  temperature  and  low  density.  Under  these  conditions, 
the  channel  conducts  nearly  all  of  the  electrical  current.  The  ambient  air  conducts 
very  little  electrical  current,  but  there  is  a  cylindrical  shock  wave  expanding  through 
it.  Flash  lamps  have  all  these  features  except  for  an  expanding  shock  wave  in  the 
air.  The  general  approach  taken  by  the  early  calculations  was  to  assume  that  the 
channel  had  a  uniform  temperature,  pressure,  current  density,  and  electrical 
conductivity.  Nearly  all  of  these  theories  assume  that  the  channel  has  an  initial 
radius  or  implies  that  there  is  one.  They  also  imply  that  the  total  mass  of  the  gas  in 
the  initial  channel  is  constant.  Various  approximations  were  then  used  to  solve  the 
equations  of  motion  of  the  gas  for  the  radius  of  the  channel  and  the  resistance  of 
the  arc  as  functions  of  time. 

Engel  et  al.1  reviewed  these  calculations  and  compared  them  with  the  experimental 
results  by  Akiyama  et  al.,2  who  measured  the  current  for  arcs  through  air  at  various 
pressures  and  arc  lengths.  The  resistances  for  the  arcs  were  then  calculated  from 
the  current  and  a  circuit  model  of  the  experiment.  The  resistance  was  high  at  the 
start  of  the  arc,  decreased  to  a  minimum  near  the  peak  current,  and  slowly  increased 
afterward.  Engel  chose  one  of  the  results  for  comparison  with  a  number  of  theories. 
Most  of  the  reviewed  theories  approximated  the  resistance  of  the  arc  up  to  about 
the  peak  current.  Some  of  the  theories  had  to  be  scaled,  but  some  of  the  scaling  may 
not  be  justified.  Most  of  the  reviewed  theories  gave  a  resistance  that  continued  to 
decrease  after  the  maximum  current.  The  few  theories  that  gave  an  increasing 
resistance  did  not  increase  as  much  as  the  experimental  result.  An  agreement  for 
the  resistance  at  all  times  was  achieved  by  combining  1)  Engel’s  proposed  equation 
for  time  dependence  of  the  radius  that  is  an  integral  of  a  function  of  the  current  and 
2)  an  equation  in  one  of  the  reviewed  theories3  that  gave  the  resistance  of  the  arc  as 
a  function  of  the  current.  A  close  examination  of  the  current  and  the  circuit  model 
used  by  Akiyama  et  al.2  to  derive  the  arc  resistance,  however,  showed  that  this  arc 
resistance  can  be  questioned. 

A  detailed  numerical  study  of  an  arc  came  from  the  research  on  thunder  and 
lightning.4  This  study  is  an  extension  of  an  earlier  finite  element  model  for  shock 
waves  from  a  line  source,5  where  the  hydrodynamic  equations  that  include  a 
pseudoviscous  term6  for  the  formation  of  a  shock  wave  were  numerically  solved. 
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These  equations  were  extended  by  adding  a  tenn  to  the  energy  balance  equation 
that  is  a  sum  of  the  energy  added  to  or  lost  from  a  finite  element4  by  various 
processes:  Joule  heating,  thermal  conductivity,  radiation  emitted  by  the  finite 
element,  and  the  radiation  absorbed  from  the  rest  of  the  finite  elements.  This  work 
was  later  reviewed  by  Akram  and  Lundgren7  and  by  a  continued  publication, 
Akram,8  which  used  a  slightly  different  set  of  finite  difference  equations  and  a 
slightly  different  equation  of  state  for  air.  The  equations  for  the  thermal 
conductivity  of  air  and  the  electrical  conductivity  of  air  are  the  same  as  Plooster’s.4 
Unfortunately,  the  equation  of  state  for  air  and  the  equation  for  the  electrical 
conductivity  of  air  in  Akram  and  Lundgren7  have  typographical  errors,  both  of 
which  are  corrected  here.  Plooster4  uses  an  equation  for  the  electrical  conductivity 
of  air  that  was  a  modification  by  Olsen9  for  atmospheric-pressure  plasmas  from  an 
expression  for  an  ionized  gas  by  Spitzer  and  Harm.10  After  referring  to  Olsen,9 
however,  the  modification  was  not  clearly  stated.  There  were  changes  in  notation 
and  possible  typographical  errors.  In  addition,  the  subject  of  radiation  transport, 
where  the  black  body  radiation  from  one  finite  element  is  absorbed  and  heats  a 
distant  finite  element,  was  not  explained  very  well  in  Plooster4  and  Akram  and 
Lundgren.7 

A  detailed  description  of  the  Lagrangian  hydrodynamic  equations  for  the  finite 
elements  are  developed  and  presented  here  and  are  different  from  those  of  Plooster4 
and  Akram  and  Lundgren.7  The  equation  of  state  for  air,  the  electrical  conductivity 
of  air,  and  a  modified  equation  for  the  pseudoviscosity  are  also  discussed  in  detail. 
It  is  assumed  that  the  energy  transport  from  one  finite  element  to  the  others  is  by 
photoionization,  where  the  ultraviolet  light  ionizes  the  neutral  atoms  in  a  finite 
element.  The  addition  of  the  energetic  electrons  and  their  density  increases  the 
electrical  conductivity  of  the  finite  element.  The  bound-bound  and  free-free 
radiative  transfer  processes  present  in  a  plasma  that  is  in  local  thermal  equilibrium 
are  ignored.11-13  It  is  assumed  that  these  processes  have  little  effect  on  the  electrical 
conductivity.  The  free-bound  process  where  a  free  electron  recombines  with  an  ion 
by  emitting  a  photon  is  included  by  using  a  rough  estimation.  These  equations  were 
solved  for  the  assumed  current  profile  used  by  Plooster4  and  Akram  and  Lundgren7 
both  when  there  is  no  photoionization  and  when  there  is  photoionization.  The 
temperature  and  pressure  at  the  very  center  of  the  arc,  the  radius  of  the  arc,  and  the 
position  of  the  shock  wave  versus  time  are  then  compared  with  Akram  and 
Lundgren’s7  results. 

2.  Hydrodynamic  Equations 

In  the  Lagrangian  description  for  a  1-D  cylindrical  coordinates,  let  the  initial  radial 
position  of  a  very  small  portion  of  the  gas  bex.  At  some  later  time,  this  small  portion 
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of  the  gas  is  moved  to  a  new  location  r(x,t).  This  portion  then  has  a  velocity  u(r(x,t)), 
a  density  p(r(x,t)),  a  pressure  p(r(x,t)),  and  an  internal  specific  energy  s(r(x,t))  in 
units  of  energy  per  unit  mass,  which  all  follow  a  set  of  differential  equations.  In 
some  of  these  equations,  however,  it  is  convenient  to  use  the  specific  volume  v  that 
is  simply  the  reciprocal  of  the  density,  v  =  \/p.  The  differential  equations  for  these 
quantities  are  found  from  the  conservation  of  mass, 


v(r(x,  t)) 


r(x,t)  dr(x,t) 

V°—  ~ 


(1) 


the  conservation  of  momentum, 


du(x,t) 

dt 


~Vn 


d(p(r(x,t))+  q(r(x,t ))  r(x,t ) 
dx  x  ’ 


(2) 


and  the  conservation  of  energy, 

Wrt,*nrlIJ"  =  -(p(p(r(*.t)),nr(*,t))  +  cKr(x,tpd-^l+W,  (3) 


where  W  is  the  power  per  unit  mass  from  Joule  heating,  radiation  loss  and  gain, 
thennal  conduction,  and  so  forth.  The  velocity  is  simply 


it(x,  t) 


dr(x,t ) 
dt  ’ 


(4) 


and  T(r(x,t))  is  the  temperature.  v0  is  the  specific  volume  of  the  ambient  air  that  is 
uniform  at  the  beginning.  The  tenn  q(r(x,t))  is  the  pseudoviscosity  term6 

q  =  p(p0cAx)2  when  ^  <  0 


or 

q  —  0  when  ^  >  0,  (5) 

where  c  is  a  constant  with  no  units  that  is  on  the  order  of  unity.  When  this  constant 
is  increased,  the  width  of  the  shock  wave  front  is  spread  out  over  a  number  of  grid 
points  and  smoothes  out  some  of  oscillations  near  the  front.  A  value  of  about  1 .4 
was  used  in  the  following  results.  This  expression  for  the  pseudoviscous  tenn  is  not 
convenient  to  use  for  the  finite  difference  equations  that  are  forward  in  time  and 
centered  in  space.14  In  this  scheme,  the  values  of  the  position  r(x,t),  velocity  u(x,t), 
specific  volume  v(x,t),  specific  internal  energy  e(x,t),  and  pressure  p(x,t)  are  used  to 
advance  their  values  to  a  short  time  step  forward  At:  r(x,t+At),  u(x,t+At),  v(x,t+At), 
e(x,t+At),  and  p(x,t+At).  The  time  derivative  of  the  specific  volume  as  it  stands 
depends  on  its  known  present  value  and  its  unknown  future  value.  It  can  be 
expressed,  however,  in  terms  of  the  known  distribution  in  space  of  variables  at  a 
given  time  by  starting  with  the  time  derivative  of  the  conservation  of  mass  equation 
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dv 

i 

dr  dr 

d2r 

1 

dr 

du 

~dt 

p0x 

_dt  dx 

+  rT^“ 
dtdx. 

p0x 

U~dx 

+  r  — 

OX 

1  d(ru ) 
Po  dx 


(6) 


to  give 

q  =  pP-il)2  (dJpl)2  when  ^<0  or  ,  =  0  when  >  0.  (7) 

“^\xj\dxj  dx  “  dx  v/ 

A  common  convention  in  writing  the  finite  difference  equations  is  to  give  a 
superscript  and  a  subscript  of  integers  to  the  variables  for  the  boundary  between  the 
elements,  as  shown  in  Fig.  1,  where  rj1  is  the  position  of  the  boundary,  uj1  is  the 
velocity  of  the  boundary,  and  fj1  is  a  flux  that  is  crossing  the  boundary.  The  flux 
can  be  a  heat  flux  or  a  radiation  flux  as  examples.  The  superscript  means  that  the 
time  for  these  variables  is  at  the  nth  time  step,  and  the  subscript  means  that  the 
variables  began  at  the  /th  starting  position.  If  the  starting  positions  are  evenly 
spaced  by  Ax,  rj1  is  equivalent  to  r(j  Ax, t„),  and  the  position  at  the  next  time  step 
ryn+1  is  equivalent  to  r(j  Ax, t„+At„),  or  r(j  Ax, tn+i)  for  the  next  time  step.  The 
starting  positions  xj,  do  not  have  a  superscript  because  they  do  not  change  with  time. 
The  positions  rj1  and  rJL+l  define  the  boundaries  for  an  element  of  gas  with  intrinsic 
properties,  which  are  conventionally  indicated  by  a  subscript  with  an  added  half 
integer  in  Fig.  1:  the  pressure  pJ+1/2,  temperature  TJn+1/2,  density  pJ+1/2,  specific 
internal  energy  sf+1/2,  and  thennal  conductivity  /r”+1/2. 


T  n 

rj-l/2'  'j-1/2; 


n  n  n 

Pj-i/2'  A  i/2*  *1-1/2 


n  T  n 

Pj+l/2'  j+1/2' 

n  n  n 

Pj+l/2'  Ej+l/2'  *j+l/2 

I 

I 

I 

I 

I 

-  6r  — *k —  6r  — 


M  -p  M 

Pj+3/2'  j+3/2' 

n  n  n 
Pj+3/2'  Ej+3/2'  Kj+3/2 


■j+1 


Xj ,  rj,  Uj 


xj+i '  rj+i'  uj+i 


Fig.  1  Notation  used  for  the  finite  difference  equations 

The  finite  difference  equations  are  used  in  the  following  sequence.  The  velocity, 
position,  and  density  are  first  stepped  forward  in  time: 

n 


n+ 1  n 

Uj  -Uj 


rj_ 

Ax  x„ 


n  n  ,  n  n 

P  j- 1/2  _  P 7+1/2  ^  j- 1/2  _  <ij+ 1/2 


71  +  1 


Tj  +  AtnUj 


71  +  1 


(8) 

(9) 
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1 


(10) 


^n+ 1 


P  y+1/2  n+ 1  P  o 


V, 


(pp-(p'f ' 


n+1  _  1  _  (*7+1)  (+/) 

^  -  v?«  -  P°  (rn S.)2_(r»«)! 


(ID 


Next,  the  pseudoviscosity  is  stepped  forward  in  time: 

v  2 


«+l 


#/+l/2 

and 


^  w+1  |  n  ^ 

P  j+m  P y+1/2 

l  2  J 

v  xj+l+xj  7 

when  <0  (12) 


=  0  when  r/^ujVV  -  r;n+1uyn+1  >  0. 


(13) 


Finally,  a  transcendental  equation  must  be  numerically  solved  for  the  new 
temperature. 


(  ^7+1/2  ’  Pj+ 1/2  )  +  P  (  Tj+ 1/2  ’  Py+1/2  ) 


/  «  +  l  «  +  l 


«+l  n+ 1  \ 


i  “O  ", 


poAx 


s{Tj+v2’Pj+m) 


n+ 1 

“#7+1/2 


f  n+ 1  «+l  «+l  «+l 

Wj+l-Vj  Uj 


(14) 


p0Ax 


KWj+m 


The  function  p(T,p)  (Eq.  30)  is  the  thermal  equation  of  state,  and  the  function  s(T,p) 
(Eq.  29)  is  the  caloric  equation  of  state  for  the  gas;  more  details  of  these  functions 
will  be  given  later.  If  the  gas  is  ideal,  these  equations  are  simple  and  easily  solved 
for  Tjl+1.  If  the  gas  is  real,  however,  these  equations  may  not  be  solved  in  closed 
form  for  the  new  temperature.  Instead,  the  old  temperature  Tj1  serves  as  a  first  guess 
for  the  new  temperature,  and  the  equation  is  tested  for  equality.  Based  on  the  results 
of  this  test,  another  temperature  is  judiciously  chosen  and  tested.  By  using  the  mean 
value  theorem  and  testing  for  equality,  an  acceptable  range  for  the  temperature  can 
be  found  in  which  there  is  equality.  A  temperature  within  this  range  is  chosen  for 
7) +1/2  to  find  the  new  specific  energy 

£7+1/2  -£[Tj+ll2’Pj+\/2  )  (13) 

and  the  new  pressure 

P 7+1/2  —  P  (  ^7+1/2  »  P 7+1/2  )  •  (16) 

These  new  values  are  then  used  for  the  next  time  step. 
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The  value  for  the  time  step  Atn  has  limits.  These  difference  equations  become 
unstable  if  Atn  is  too  large.  If  At,,  is  too  small,  the  computation  time  for  a  total 
solution  over  a  range  of  time  increases.  Plooster15  used  the  stability  relation  given 
by  Richtmyer  and  Morton,16 


At,-  <  —  ^ 


J  4 a2  Uj-Uj+1  ’ 


(17) 


where  a  is  a  parameter  that  originated  from  the  pseudoviscosity.  This  condition  has 
a  simple  interpretation.  The  numerator  is  the  spacing  between  2  positions.  The 
denominator  is  the  closing  velocity  of  the  2  positions  when  uj  >  Uj+i.  Thus,  this 
fraction  is  the  additional  time  when  the  distance  between  the  2  positions  is  zero  if 
they  continue  to  move  at  their  velocities.  This  cannot  happen  because  the  positions 
represent  material  that  cannot  pass  through  another  and  result  in  a  negative  width 
for  the  element.  If  both  velocities  are  zero  or  if  Uj<  Uj+i  for  an  expanding  element, 
then  the  time  step  is  taken  to  be 


At:  < 


4a2 


5  7+1/2 


(18) 


where  Sj+1j 2  =  yjYPj+1/2  /  d;+ 1/2  is  the  local  speed  of  sound  with  y  =  1.4.  Because 
each  element  has  a  different  time,  all  the  elements  must  be  examined  for  a  minimum 
time,  and  the  time  step  is  set  to  a  fraction  of  this  time.  A  value  for  a  that  is  taken  to 
be  independent  here  from  the  constant  c  in  Eq.  5  can  be  found  with  some  trial  and 
error,  but  a  =  1  was  sufficient  for  the  following  calculations. 


In  the  conservation  of  energy  equation,  the  W  term  is  the  sum  of  the  energy  per  unit 
time  per  mass  that  is  added  or  lost  by  the  element  Wj+1/2  =  W joule  ~  Wbb  + 
W therm-  W joule  is  the  Joule  or  ohmic  heating. 


WJoule  =  E 


2  <Jj+ 1/2 
P  j+V2 


(19) 


E  is  the  common  electric  field  for  all  the  elements  that  is  given  by  the  total  current 
divided  by  the  sum  of  the  products  of  the  electrical  conductivity  and  the  cross 
sectional  area  of  all  the  elements. 


E  = 


/(t) 


(20) 


where  a,-  =  7r( rf+1  —  rf  ')  and  the  total  current4  7  is  assumed  to  be 

/(t)  =  4.0  x  104(exp(— 4.0  x  104t)  —  exp(— 4.0  x  105t)). 


(21) 
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The  maximum  for  this  current  is  about  27.9  kA  at  about  6.4  (.is.  The  power  loss 
from  the  black  body  radiation  according  to  Plooster4  is  Qbb  —  4asbT4kb.  The  tenn 
4 osbTA  is  the  integral  of  the  spectral  radiation  intensity  over  all  the  frequency  and 
solid  angle,  where  asb  is  the  Stefan-Boltzmann  constant  and  kb  is  an  absorption 
constant  that  was  arbitrarily  chosen  to  be  1 .0  cm"1  in  the  gas  at  the  ambient  density. 
Thus,  the  specific  power  loss  per  unit  mass 


Wbb  = 


4  osbT*  1.0 
Po 


(22) 


Wthemi  is  the  heat  conducted  into  an  element  from  its  neighbors.  According  to 
Akram  and  Lundgren,7  W therm  in  differential  fonn  is 


W therm 


(23) 


k  is  the  heat  conductivity  of  the  element.  Casting  this  differential  expression  into  a 
finite  difference  form  is  not  obvious  because  it  calls  for  a  differentiation  of  a 
product  of  the  intrinsic  quantities  T,  p,  and  k  with  a  spatial  variable  x.  Transforming 
the  partial  derivatives  with  respect  to  x  into  the  partial  derivatives  with  respect  to  r, 
however,  gives  a  simpler  differential  that  is  easier  to  interpret.  Using  the 
conservation  of  mass,  the  operator 


5  _  p0  x  d 
dx  p  r  dr 


(24) 


is  used  to  express  Eq.  23  as  the  familiar  radial  component  of  the  heat  conduction 
equation  in  cylindrical  coordinates 


W, 


therm 


15/  5T\1 

rdr  \Kr  dr)\’ 


(25) 


where  the  expression  in  the  parenthesis  is  the  radial  component  of  the  heat  flux 
vector  h  —  —kVT,  and  the  expression  in  the  brackets  is  the  radial  component  of  the 
divergence  of  the  heat  flux.  This  suggests  a  way  to  translate  this  equation  to  the 
finite  elements.  Let  hj  =  k  J~1/2Sr  7+1/2  be  the  heat  flux  entering  the  element  on  the 

left,  where  it  is  assumed  that  the  temperature  at  rj  is  Tj.in  and  the  temperature  at  rm, 
rm  =  (rj+ 1  +  rj)/ 2,  is  T,  m.  The  heat  flux  leaving  the  element  on  the  right  is  hj+1  = 

K  ‘ 1+1/2 ^ 1+3/2 ,  where  it  is  assumed  that  the  temperature  at  rj+ 1  is  7)+ 3/2.  After 
multiplying  the  fluxes  by  their  areas,  the  total  heat  entering  the  element  is 


Q  =—r 

■ J  Sr,  L 


2nrj{Tj-\n 


-Tj+m)-2^rj+x{Ti+v2 


(26) 


This  is  then  divided  by  the  mass  within  the  element  for  the  specific  heat  Wtherm, 
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w. 


therm 


pj2nrm28r  pj4nrmSr2 


[2nrj(Tj_1/2  7}+ 1/2)  2n rj+1(Tj+1/2  7j+3/2)].  (27) 


After  applying  77  =  rm  -  Sr  and  77+1  =  rm  +  Sr,  this  can  be  written  as 


W therm 


2  P  / 


T.  3+7\  1—2 7\  1  T.  3 ~T .  i 
'  '  2  ;  2  rt  \5  _|_  ^+2  f  2 


5r2 


rmSr 


(28) 


The  thennal  conductivity  in  units  of  erg  cm"1  sec'1  deg'1  is  a  sum  of  3  terms. 

dA 

k  =  0.1488crr  +  2.2xl013p-^  +  1000(l-241)r1/2 


4,7 


(29) 


The  first  tenn  is  the  thermal  conductivity  by  the  electrons  detennined  from  the 
Wiedermann-Franz  law,  where  a  is  the  electrical  conductivity.  The  second  term  is 
from  the  diffusion  of  the  recombination  energy,  where  Ao  is  from  Eq.  32 
(Section  3).  The  last  term  is  the  contribution  from  the  neutral  particles,  where 
(1  ~Ai)  is  the  fraction  of  neutral  particles  and  Ai  is  also  from  Eq.  32  (Section  3). 

The  correct  initial  conditions  for  an  arc  at  the  very  beginning  of  a  current  pulse  are 
speculative,  but  Plooster4  proposed  2  different  initial  conditions.  Even  though  these 
conditions  may  not  be  correct,  it  is  hoped  that  the  results  would  not  depend  on  the 
exact  initial  conditions  at  later  times.  One  initial  condition  has  a  channel  of  gas  at 
ambient  pressure  with  an  assumed  radius  and  high  temperature.  These  conditions 
may  occur  just  prior  to  the  main  discharge  of  a  lightning  stroke  when  a  large  current 
flows  between  the  ground  and  the  cloud.  Just  before  the  main  discharge,  there  are 
leading  strokes  emitting  from  the  ground  to  the  cloud  and  other  leading  strokes 
emitting  from  the  cloud  to  the  ground.  These  leading  strokes  propagate  toward  each 
other  until  they  make  contact,  initiating  the  main  discharge.  This  also  occurs 
between  clouds.  It  is  assumed  that  these  leading  strokes  heat  a  channel  of  air  at 
atmospheric  pressure.  The  other  initial  condition  has  a  channel  of  gas  at  the  ambient 
density  with  an  assumed  radius  and  high  temperature,  which  may  occur  when  an 
arc  is  formed  so  quickly  that  the  gas  has  yet  to  move.  This  condition  implies  that 
the  gas  was  heated  by  a  current  for  some  time  before  the  gas  is  allowed  to  move. 
Plooster4  and  Akram  and  Lundgren7  ignore  this  initial  current,  and  simply  start  from 
zero  current  and  allow  the  gas  to  move  at  the  same  time.  From  Akram  and 
Lundgren’ s  graphs,7  the  initial  temperature  for  the  channel  was  15,000  K  but  its 
initial  radius  was  not  stated.  An  initial  radius  of  0. 1  cm,  however,  gave  comparable 
results. 

Using  Akram’s8  and  Plooster’s4  equations  for  the  electrical  conductivity  of  air  (Eqs. 
36  and  37;  Section  3)  Akram’s8  corrected  equation  of  state  for  air  (Eqs.  30,  3 1,  and 
32;  Section  3),  and  assuming  that  the  channel  is  air  at  ambient  density  uniformly 
heated  to  a  temperature  of  15,000  K  within  a  radius  of  0.1cm,  these  finite 
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difference  equations  were  solved  without  radiation  transfer.  At  the  time  6.4  ps  for 
the  maximum  current  27.9  kA,  the  results  gave  the  distribution  of  the  temperature 
(Fig.  2),  pressure  (Fig.  3),  density  (Fig.  4),  and  current  density  (Fig.  5).  It  is 
remarkable  that  all  of  these  quantities  are  fairly  unifonn  within  the  channel  radius 
of  about  0.82  cm  and  that  a  shock  wave  that  is  the  peak  in  the  pressure  at  1 . 1  cm  in 
Fig.  3  is  already  at  some  distance  from  the  channel.  This  supports  the  assumptions 
used  in  the  earlier  models. 


Position  (cm) 


Fig.  2  Temperature  vs.  position  at  6.4  ps 


Fig.  3  Pressure  vs.  position  at  6.4  ps 
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Fig.  4  Density  vs.  position  at  6.4  ps 


16 

14 

12 

T  10 

1 

f 

1 

>■ 

•|  8 

Q 

1 

3  6 

1 

4 

1 

,  V 

Position  (cm) 

Fig.  5  Current  density  vs.  position  at  6.4  ps 

3.  Equation  of  State  for  Air 

The  pressure  of  air  as  a  function  of  its  temperature  and  its  density  is  the  thermal 
equation  of  state,5 

p(T,p)=pRT[l+A0  +  2(A1+A2)],  (30) 

and  the  specific  internal  energy  of  air  in  ergs  per  gram  is  calculated  by  the  caloric 
equation 

£(T,p)  =  RT(S  +  A0/2  +  3(A1+A2-))+A0I0  +  A1I1+A2I2.  (31) 
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In  these  equations,  An  is  the  fraction  of  the  molecules  that  are  dissociated:  A o  =  0 
when  there  are  no  dissociated  molecules  and  Ao=  1  when  all  the  molecules  are 
dissociated.  Ai  is  the  fraction  of  the  atoms  that  are  singly  ionized:  A  /  =  0  when  none 
of  the  atoms  are  singly  ionized  and  Aj  =  1  when  all  the  atoms  are  singly  ionized. 
A 2  is  the  fraction  of  the  singly  charged  ions  that  become  doubly  charged:  A 2  =  0 
when  none  of  the  singly  charge  ions  is  doubly  charged  and  A 2  =  1  when  all  of  the 
singly  charged  ions  are  now  doubly  charged.  Ao,  Ai,  and  A 2  are  functions  of  the 
temperature  and  density  found  by  applying  the  Saha  equation13  for  the  local  thermal 
equilibrium  between  2  states  separated  by  an  energy.  In  this  case,  the  lower  energy 
state  for  Ao  is  the  molecule,  and  the  higher  energy  state  is  when  the  molecule 
dissociates  into  atoms.  The  lower  energy  state  for  A 1  is  the  atom,  and  the  higher 
energy  state  is  the  singly  charged  ion.  The  lower  energy  state  for  A2  is  the  singly 
charged  ion,  and  the  higher  energy  state  is  the  doubly  charged  ion.  After  correcting 
a  typographical  error,  these  fractions,  according  to  Akram  and  Lundgren,7  are 


A 

A 

A 


=  2 


l  +  yjl  +  2BQ 


Co  p 


exp\ 


0 


At  l-exp (-T/T) 


=  2 


l  +  Jl  +  2  Bl 


7-3/2 


exp 


'jC 

2  WlT 


=  2 


1  +  B„  +  Jl  +  +  Bl 


_  2 


C,P 


t3/2 


exp 


291 T 


(32) 


The  values  Io,  Ii,  h  and  Co,  Ci,  C2  in  the  Table  are  from  Akram8  for  an  air  that  has 
properties  that  are  a  weighted  sum  of  the  corresponding  property  of  a  gas  with  a 
ratio  in  number  of  0.788  N2  and  0.212  O2.  Compounds  of  nitrogen  and  oxygen  and 
their  ions  that  may  form  in  and  around  the  arc  are  ignored:  NO,  NO+,  O3,  and  so 
forth.  The  trace  atmospheric  gasses  are  also  ignored:  argon,  carbon  dioxide,  and  so 
forth. 


1 
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Table  Constants  for  the  equation  of  state  from  Akram.8  The  gas  constant  per  gram  is  the 
universal  gas  constant  divided  by  the  molecular  weight. 


Parameter  name 

Symbol 

Value 

Units 

Molecular  weight 

Mw 

28.86 

g/mole 

Gas  constant  per  gram 

R 

2.880946  x  106 

erg/K/g 

Specific  dissociation  energy 

10 

2.930557  x  1011 

erg/g 

Specific  first-ionization  energy 

11 

9.591766  x  1011 

erg/g 

Specific  second-ionization  energy 

12 

2.062056  x  1012 

erg/g 

Saha’s  constant  for  dissociation 

CO 

4.817588  x  10-1 

cm3K3/2/g 

Saha’s  constant  for  first-ionization 

Cl 

1.101416  x  107 

cm3K3/2/g 

Saha’s  constant  for  second-ionization 

C2 

1.033965xl07 

cm3K3/2/g 

Vibration  temperature 

Tv 

3150 

K 

Ambient  temperature 

To 

300 

K 

Ambient  pressure 

po 

1.01325  x  106 

dyne/cm2 

Ambient  density 

po 

1.172358  x  IQ-3 

g/cm3 

The  electrical  conductivity  of  the  gas  depends  on  the  number  density  of  the 
electrons  ne,  number  density  of  the  neutral  “air  atoms”  no,  number  density  of  the 
singly  charged  ions  m,  and  number  density  of  the  doubly  charged  ions  H2.  These 
number  densities  are  related  to  the  fractions  Ao,  A 1,  and  A 2  but  they  are  not  for  a 
given  number  density  of  “air  atoms”  N  where  N  =  Av  p/14.43.  A  v  is  Avogadro’s 
number,  and  14.43  is  half  of  the  molecular  weight  for  the  “diatomic  gas  Ain”: 
28.86.  In  terms  of  these  fractions,  the  number  densities  are  nm  =  (1  -  Ao)N/2  for 
molecules,  na  =  Ao(\-Ai)N  for  atoms,  no  =  2nm  +  na  =  (1  -AoAi)N  for  the  total 
number  of  atoms  that  are  associated  in  molecules  and  atoms  that  are  dissociated, 
m  =AoAi  (1  —Ai)  N  for  singly  charged  ions,  and  m=  A0A1A2N  for  doubly 
charged  ions.  Note  that  N  =  no  +  m  +  m.  Thus,  the  number  of  electrons  is 
ne  =  1  ni  +  2  H2  =AoAi  (1  +  4y)7V,  and  the  total  number  of  ions  is 
n+  =  m  +  ni  =  Ao  A iN.  The  mean  ionic  charge  that  will  be  used  in  the  equation  for 
the  electrical  conductivity  is  defined  as 


Z  = 


L2>  A . 


ne~ 


(33) 


where  Zj  is  the  charge  of  the  ion  ( Zi=  1,  Z?  =  2,  etc.)  and  n/s  are  their  respective 
number  densities.  In  tenns  of  the  fractions,  the  mean  ionic  charge  up  to  the  doubly 
charged  ions  is 


Z 


1 «,  +  2 1  n2 

ne 


A0A1(l-A2)  +  4A0AlA2  (1  +  3  A2) 

A0Al(l  +  A2)  (1 +  A2) 


(34) 
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Because  of  the  energy  dependence  of  Ao,  Ai,  and  A 2,  most  of  the  atoms  are  singly 
charged  before  there  are  doubly  charged  ions.  Also,  most  of  the  molecules  are 
dissociated  into  atoms  before  there  are  singly  charged  ions.  Thus,  the  following  are 
very  good  approximations  that  are  used  by  Plooster4  and  Akram  and  Lundgren7: 
ne  =  A0Ai(l  +  A2)N ~  (Al  +  A2)N,  n+  =  A0AiN~  A,N,  no  =  (1  -A0A,)N~  (1  -Ai)N, 
and 


(l  +  3  A2)  Al  +  3A2 
(l  +  A2)  At  +  A2 


(35) 


4.  Electrical  Conductivity  of  Gasses 


The  electrical  conductivity  used  by  Plooster4  and  Akram  and  Lundgren7  in  Siemens 
per  centimeter  (S/cm)  is 


4.173xl(Tl0(4  +  4)r~2 
2.0xl0-|5(l-4)  +  Jl(ai))/ 


(36) 


where  the  average  electron-ion  cross  section  ( at)Av  is  given  by 


(fl/L  = 


2.8xl(T 


+  3A2 
Aj  +  A2 


Ln 


1.727  xl0“5r 


(Ap) 


1/3 


A\  A2 
A ^  +  3A2 


(37) 


Using  the  approximations  for  ne,  no,  and  n+  from  the  equation  of  state,  these 
equations  are 


4.173xl(T10ng/2 
2.0xl0“15n  +n  ,(a.)A 

n  +  \  (  /  Av 


(38) 


and 


(*/L 


2.8  x  10-6Z2 

F2 


Ln 


1.727  x  10_5T  ] 

(ap)i/3z  r 


(39) 


Plooster4  simply  states  that  this  equation  is  a  modification  by  Olsen9  of  the  electrical 
conductivity  for  an  ionized  gas  of  Spitzer  and  Harm10  for  atmospheric-pressure 
plasmas.  Olsen  did  not  provide  any  explanation  for  the  modifications.  Also,  there 
were  some  typographical  errors,  changes  in  notation,  choice  for  some  terms,  and 
possible  omissions.  Thus,  another  equation  for  the  electrical  conductivity  will  be 
developed  here. 


The  electrical  conductivity  of  gasses10,17  in  electrostatic  units  is 


a  — 


e2ne 


mCZjnjQj 
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(40) 
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To  convert  the  unit  of  this  electrical  conductivity  to  Siemens  per  centimeter,  divide 
it  by  c2xl0"9,  where  c  is  the  speed  of  light  in  centimeters  per  second  or  by 
multiplying  it  by  47is0.  e  is  the  charge  of  the  electron  in  Statcoulombs,  ne  is  the 
number  of  free  electrons  per  cubic  centimeter,  m  is  mass  of  the  electron  in  grams, 
and  C  is  the  root  mean  square  of  the  electron  speed  having  a  Maxwell-Boltzmann 
speed  distribution 


C  = 


(41) 


where  k  is  the  Boltzmann  constant  in  ergs  per  Kelvin.  Substituting  this  into  the 
Eq.  4 1  gives 

e2  nP  „„  nP 

a  = - - - t — - =  4.179  x  10“10  -t — -  (42) 

c2  x  1 0  _9V3mk  T7.YijnjQj  T^njQj 

for  the  conductivity  in  Siemens  per  centimeter.  The  summation  is  over  the  products 
of  the  number  density  for  a  component  of  the  gas  and  the  cross  section  for  electron 
collision  with  the  component.  The  components  for  air  would  be  molecules,  atoms, 
and  ions  of  nitrogen,  oxygen,  and  so  forth.  The  summation  is  broken  into  2  terms. 

YnjQJ=2.0xlOl5n0  +  YnJQJ.  (43) 

7=0  7=1 


The  first  tenn  is  the  assumed  cross  section  in  square  centimeters  for  the  electrons 
colliding  with  neutral  atom  or  molecule  multiplied  by  the  number  of  “air  atoms” 
even  when  they  are  bonded  in  molecules.  This  term  implies  that  the  cross  section 
for  a  diatomic  air  molecule  is  about  twice  that  of  the  atom  in  the  molecule,  which 
is  a  reasonable  approximation.  The  second  term  is  the  product  of  the  cross  section 
and  the  number  density  of  the  ions  having  a  charge  j:  j  =  1  for  singly  charged  ions, 
y  =  2  for  doubly  charged  ions,  and  so  forth.  The  second  tenn  is  found  by  assuming 
that  all  of  the  atoms  are  completely  ionized  and  equating  the  conductivity  to 
electrical  conductivity  of  an  ionized  gas  according  to  Spitzer  and  Hann10: 


cr  = 


2  mC3 


e1  Z\n[qC2)V?>7r 


,3/2 


r*(z), 


(44) 


where  Z  is  the  mean  ionic  charge  (Eq.  31).  After  equating  the  2  conductivities, 
Eqs.  40  and  44  with  no  =  0,  the  summation  is 


Yrfij 


eAne  Z  \n(qC2) 


2/^*4 


2m  C 


3/r 


v  z.  j 
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rJz)' 


(45) 
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After  substituting  the  expression  for  the  root  mean  square  of  the  velocity  and 
collecting  the  constants, 


Yjifi 


(  4 

e  \ 

f  3^7 

3/2^ 

1  ne  Zln (<?C2) 

18k2! 

V 

9 

rE(z)  t 1 

and,  evaluating  them, 


1.58689  x  IQ6 

rE{z) 


ne  Zln ) 

rjn2 


(46) 


(47) 


The  yE  (Z)  comes  from  the  electron-electron  collisions,  and  it  too  is  a  function  of 
the  mean  ionic  charge  Z18,19: 


rE(z) 


"3;r)l  +  2.966Z  +  0.753Z2 
v32~ J  1  +  1.198Z  +  0.222Z2  ' 


(48) 


A  ye  =  0.567  would  give  the  2.8  x  10'6  coefficient  in  Eq.  37,  which  is  close  to  0.574 
when  Z  =  1  in  Eq.  48.  The  value  for  qC2  depends  on  the  integration  over  the  impact 
parameter  of  the  collision.  The  upper  limit  of  the  integral,  however,  logarithmically 
approaches  infinity  and  must  be  cut  off  at  some  distance.  The  lower  limit  of  the 
integral  is  undefined  for  a  zero-impact  parameter.  Cohen  et  al.20  used  the  Debye’s 
shielding  radius  for  the  upper  limit  where  the  distribution  of  the  electrons  and  ions 
beyond  this  distance  has  little  influence  on  the  collision.  The  lower  limit  of  the 
integral  is  taken  to  be  the  classical  distance  for  closest  approach  of  the  electron  to 
the  ion.  Evaluating  the  integral  with  these  limits  gives 


qC2  = 


3 


(kT\ 


3/2 


(49) 


Thus,  the  electrical  conductivity  is 

_ _ 4.179  x  icr10ngr3/2 _ 

(2.0  x  10  15n0Z2  +  1.158689  x  \0~6neZln(qC2)/ yE  (Z))  ’ 


(50) 


or  when  using  the  fractions  from  the  equation  of  state, 

= _ 4.179xlQ-10Z0Z1(l  +  Z2)r3/2 _ 

^  (2.0xl0'15r2(l-fioZl)  +  1.158689xl0“6Z0Z1(l  +  3Z2)/«(^C2)/yi?(Z))  ’ 

where  qC2  and  Ye(Z)  are  given  in  the  Eqs.  48  and  49. 

If  the  interionic  distance  is  used  as  the  cut-off  distance,  instead  of  the  Debye 
shielding  radius,  the  term  would  be 
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for  singly  charged  ions,  as  cited  by  Cohen  et  al.20  and  used  in  Eqs.  37  and  39.  The 
average  cross  section  for  Plooster’s4  Eq.  39  has  the  square  of  the  mean  ionic  charge 
Z,  but  Eq.  46  has  the  mean  ionic  charge  to  the  first  power.  Since  the  mean  ionic 
charge  is  close  to  1 .0  for  most  of  the  conditions  considered  here,  this  difference  is 
inconsequential.  Figure  6  compares  the  electrical  conductivities  of  the  2  equations 
as  function  of  the  temperature  for  various  densities  that  occur  in  the  calculations 
O.Olpo,  1.0po,  and  10.0po. 


Temperature  (K) 


Fig.  6  Electrical  conductivity  by  Eq.  36  ( - )  and  Eq.  50  ( - ) 

5.  Photoionization 


The  photoionization  of  air  by  a  black  body  radiator  is  presented  by  a  simple 
example.  Let  a  black  body  radiator  be  a  cylinder  with  a  temperature  Tb  and  a  radius 
Rb  that  has  a  layer  of  air  next  to  it.  The  layer  of  air  has  a  temperature  Ta,  density  p, 
and  thickness  s.  Assume  that  the  temperature  of  the  black  body  and  the  air  are 
initially  the  same.  At  some  time  let  the  temperature  of  the  black  body  be  suddenly 
increased  to  a  higher  temperature  when  the  2  regions  are  no  longer  in  thermal 
equilibrium.  The  photons  from  the  black  body  can  now  ionize  more  air  to  increase 
the  number  density  of  free  electrons,  and  these  electrons  or  photoelectrons  have  an 
energy  distribution.  Thus,  the  photoelectrons  can  change  the  electrical  conductivity 
of  the  air,  and  this  change  depends  on  time.  The  production  rate  of  the 
photoelectrons  and  their  energy  distribution  can  be  estimated  by  looking  at  the 
details  for  photoionization. 
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When  a  photon  that  has  energy  (hv)  greater  than  the  ionization  potential  of  an  atom 
(ip)  is  absorbed  by  the  atom,  the  atom  ejects  a  photoelectron  with  energy  (hv-ip)  to 
become  a  singly  charged  ion  in  its  ground  state,  h  is  Planck’s  constant,  and  v  is  the 
frequency  of  the  photon,  P(hv)  +  A  —>  A'  +  e(hv_ip) . 

The  probability  of  a  photon  being  absorbed  by  an  atom  or  molecule  and  producing 
a  photoelectron  is  related  to  a  cross  section  that  depends  on  the  energy  of  the 
incident  photon.  It  is  possible  that  the  singly  charged  ion  may  be  in  an  excited  state 
after  ejection  of  the  photoelectron  with  even  less  energy.  It  is  also  possible  that  an 
atom  absorbs  a  photon  at  some  energy  to  become  an  atom  in  an  excited  state  that 
later  decays  by  ejecting  an  electron  (autoionization).  For  air,  however,  a  photon 
may  also  be  absorbed  by  a  molecule  that  has  additional  mechanisms  to  produce  a 
photoelectron  with  less  energy.  The  molecule  may  become  a  molecular  ion, 
dissociate  into  an  atom  in  the  ground  state  or  an  excited  state  and  an  ion,  and  so 
forth.  These  latter  mechanisms  are  not  specifically  addressed  but  they  are  indirectly 
included  here  by  using  the  total  cross  section,  which  includes  all  these  mechanisms. 
Unfortunately,  the  literature  gives  these  cross  sections  in  different  units  as  functions 
of  the  frequency,  or  the  energy,  or  the  wave  number,  or  the  wave  length  of  the 
photon.  Also,  the  cross  sections  are  given  as  smooth  line  graphs,  data  point  graphs, 
or  listed  in  tables. 

In  the  following  sections,  the  cross  sections  for  the  photoionization  of  the  nitrogen 
molecule  (N2),  nitrogen  atom  (N),  oxygen  molecule  (O2),  and  oxygen  atom  (O)  are 
presented  in  units  of  Megabarn  (Mb)  as  functions  of  the  photon  energy  in  electron 
volts  (eVs).  One  Mb  is  10' 18  cm2.  Linear  interpolation  was  used  to  find  the  cross 
sections  on  a  common  energy  scale,  as  shown  in  Fig.  7.  The  total  photoionization 
cross  section  for  N221  and  O222  is  the  production  of  a  photoelectron  with  any  energy 
leaving  the  molecule  in  any  possible  state:  molecular  ion  in  any  of  its  states, 
dissociated  atom  and  ion,  and  so  forth.  The  total  photoionization  cross  section  for 
N23  and  O24  includes  the  ions  in  an  excited  state.  Some  structure  in  these  cross 
sections  due  to  the  autoionization  of  a  highly  excited  atom  was  ignored.  Note  that 
the  cross  sections  for  the  molecules  are  approximately  twice  that  of  their 
corresponding  atomic  cross  sections. 
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Fig.  7  Total  photoionization  cross  sections,  N2  ( - ),  O2  ( - ),  N  ( - ),  and  O  ( - ) 

These  cross  sections  were  used  to  find  the  energy  distribution  of  the  photoelectrons 
when  the  source  of  the  photons  is  a  black  body  radiator.  This  energy  distribution 
may  be  related  to  a  temperature.  As  an  example  in  finding  this  energy  distribution, 
assume  that  the  energy  of  the  photons  is  25  eV,  and  assume  that  the  air  is  all 
dinitrogen  and  dioxygen.  If  the  photoelectron  came  from  dinitrogen,  its  energy 
would  be  25  eV  minus  the  ionization  potential  of  dinitrogen  (15.60  eV),  or  9.40  eV. 
If  the  photoelectron  came  from  dioxygen,  however,  its  energy  would  be  25  eV 
minus  12.07  eV,  the  ionization  potential  of  dioxygen,  which  would  be  12.93  eV. 
There  will  be  photoelectrons  with  less  energy  if  the  molecular  ion  is  excited  or  if 
the  molecular  ion  dissociates.  These  photoelectrons  are  ignored,  since  the  number 
of  these  lesser  energetic  photoelectrons  are  smaller,  because  the  cross  sections  are 
smaller.  The  energies  are  weighted  by  their  respective  cross  section  and  by  their 
fractional  number  in  air  to  find  their  average  energy. 


0.788 (hv  - 1 5.60)  QN2  (, hv )  +  0.212  (hv  - 12.07)  Q02  ( hv ) 

Q(hv) 


,  (53) 


where 


Q(hv)  =  0.788  Qm(hv)  +  0.212  Qm(hv ) 


(54) 


is  an  average  cross  section.  The  0.788  factor  is  the  numerical  fraction  for  the 
dinitrogen  having  a  cross  section  Qn2(Iiv),  and  0.212  is  the  numerical  fraction  for 
the  dioxygen  having  a  cross  section  Qo2(hv).  The  trace  gasses  in  air,  argon,  carbon 
dioxide,  and  so  forth,  are  excluded:  If  the  photon  energy  is  less  than  an  ionizing 
potential,  the  respective  cross  section  is  zero. 
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When  the  air  is  exposed  to  the  black  body  radiator,  the  production  rate  of 
photoelectrons  with  this  average  energy  within  a  small  range  of  energies  d(hv)  per 
unit  volume  per  unit  time  is 

E(hv)  =  Eavg  ( hv)NmQ(hv)P(Tb,hv)d  (, hv ) ,  (55) 

where  Nm  is  the  number  density  of  the  dinitrogen  and  dioxygen  of  this  reduced  air. 
At  a  temperature  of  300  K  and  a  density  1 .  1 768 1  x  1  O'1  g/cm3,  the  number  density 
Nm,  is  2.44629x  1 019  1/cm3.  The  total  number  of  photons  within  the  energy  range 
d(hv)  that  are  emitted  from  a  unit  area  per  solid  angle  from  the  black  body’s  surface 
at  a  temperature  of  Tt  per  unit  time  is 


hv 


exp 


dv 

r  hv  ' 

,kTb; 


1 


(56) 


The  expression  in  the  brackets  is  for  the  energy  emitted  in  a  frequency  interval  dv 
per  unit  area  per  unit  time.  Dividing  this  by  the  energy  of  the  photon  hv  gives  the 
number  of  photons  per  unit  area  per  unit  time,  which  can  then  be  expressed  as  a 
function  of  the  photon  energy. 


P{hv,Tb) 


2 1.3 


c~h 


(hv)2 


exp 


rkv_' 

\kTtj 


-1 


d  (hv) 


Thus  the  overall  average  energy  is 


(57) 


E(Th) 


J hvEavg(hV)Q(ltV)  PihV’Tb)  CP,lV) 

{ hvQ(hv)  p(hv,Tb)  d {hv) 


(58) 


The  Nm  factor  cancels  out  in  this  ratio.  This  energy  average  is  then  equated  to  a 
temperature  by 


E(T„) 


m  v2  m  f  2>kTe  N 

2  2  \  m  j 


(59) 


for  the  equivalent  temperature  for  the  photoelectrons  even  though  their  energy 
distribution  is  not  a  Maxwell-Boltzmann  distribution. 


2E(T„) 

3k 


(60) 
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The  integrals  used  the  published  cross  sections  (Fig.  7),  up  to  a  photon  energy  of 
50  eV,  where  the  photon  number  flux  is  small  even  when  the  black  body 
temperature  is  60,000  K.  The  results  are  shown  in  Fig.  8  along  with  the  result, 
where  it  is  assumed  that  the  air  is  all  atomic.  For  all  atomic  air,  the  ionization 
potential  for  N  (14.53  eV),  O  (13.62  eV),  and  their  total  photoionization  cross 
sections  were  used  to  find  the  equivalent  temperature  for  their  photoelectrons.  The 
black  body  temperature  is  included  for  comparison. 


Black  Body  Temperature  (K) 


Fig.  8  Molecular  air  ( - ■),  atomic  air  ( - ),  and  black  body  temperature  ( - ) 

Since  electron  temperature  is  nearly  equal  to  the  black  body  temperature  when  the 
air  is  all  atomic  or  all  molecular,  the  electron  temperature  will  be  set  equal  to  the 
black  body  temperature  for  the  electrical  conductivity. 

The  production  rate  of  the  free  charges  in  the  thick  layer  is  found  by  first  defining 
a  small  volume  that  is  a  right  cylinder  with  a  cross  sectional  area  a  and  a  length  d. 
At  one  end  of  the  cylinder,  let  a  number  of  photons  Np(0,hv)  with  an  energy  hv 
per  energy  bandwidth  d(hv),  per  unit  area,  and  per  unit  time  cross  perpendicular  to 
the  area  a.  If  the  layer  has  a  uniform  number  density  of  neutral  atoms  or  molecules 
N,  the  production  rate  for  the  photoelectrons  per  photon  energy  in  this  volume  is 

Ne(hv)=Np(0,hv)NQ(hv)aS ,  (61) 

where  Q  is  a  cross  section.  Since  for  every  photon  that  is  absorbed  there  is  a 
photoelectron,  the  number  of  photons  leaving  the  cylinder  at  the  other  end  is 
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for 


Np(S,hv)a  =  Np(0,hv)a-Np(0,hv)NQ(hv)a  S 


(62) 


Np(S,hv)-Np(0,hv ) 
£ 


Np(0,hv )  N  Q(hv). 


(63) 


Taking  the  limit  of  S  approaching  zero  on  the  left  side,  this  equation  becomes 

dNn(hv)  .  .  .  .  .  N  ( hv )  ,  . 

- ,^rL  =  Np(0 ,hv)  N  Q(hv)=  -^-±  =  hXhv)  (64) 

do  ad 

for  the  number  of  photo  electrons  being  produced  per  unit  volume  per  unit  time. 
Finally,  let  the  number  distribution  of  the  photons  be  proportional  to  the  number 
distribution  of  the  black  body  photons  that  can  photoionize. 

Np(0,hv)  =  a  P(hv,Tb), 


where  a  includes  an  attenuation  factor  and  a  solid  angle  factor.  Thus,  a  production 
rate  factor  can  now  be  defined  as 


F(Tb)=\Q(hv)P(hv,Tl)d(hv )  (65) 

ln’o 


and 

Ae(Tb)  =  NaF(Tl)  (66) 

for  the  total  photoelectron  production  rate  per  volume,  where  N  is  the  number 
density  of  the  neutral  atoms  even  when  they  are  molecules.  Figure  9  shows  that 
F(Tb )  when  it  is  assumed  that  the  air  is  all  molecular  and  when  it  is  assumed  that 
the  air  is  all  atomic.  In  both  cases,  F(Tb )  is  almost  zero  when  the  black  body 
temperature  is  less  than  about  15,000  K  and  nearly  equal  at  higher  temperatures. 
Since  it  is  expected  that  most  of  the  photoionization  will  occur  when  the  air  is 
mostly  molecular,  however,  the  production  rate  for  the  molecular  air  was  chosen 
and  fitted  to  a  polynomial  for  temperatures  greater  than  15,000  K, 

F(x )  =  -4.754  x  10~16x5  +  5.325  x  10'ux4  -  9.127  x  10'V  +  1.664  x  10“V  +  2.774x ,  (65) 

where  x  =  Tb-15,000  when  Tb  is  greater  than  15,000  K.  F(Tb )  =  0  when  Tb  is  less 
than  15,000  K. 
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Fig.  9  Production  rate  for  molecular  air  ( - )  and  atomic  air  ( - ) 


Now  assume  that  the  thickness  of  a  cylindrical  element  in  the  arc  and  the  mean  free 
path  of  the  photons  are  small  compared  with  the  radius  of  the  black  body  Rb  so  that 
the  cylindrical  surface  of  the  black  body  can  now  be  taken  to  be  a  planar  surface. 
Let  the  surface  of  a  black  body  have  an  area  dAb  where  the  photons  are  being 
emitted  with  a  cos  6  angular  distribution  and  9  is  the  angle  between  the  normal  of 
the  area  dAb  and  the  travel  direction  of  the  photons.  The  total  production  rate  for 
the  photo  electrons  in  a  small  volume  located  at  a  perpendicular  distance  s  from  the 
black  body  surface  is 


he{Tb,s)  =  NF(Tb) 


R. 


v  Rb  +  5  j 


exp  (-kr)  cos  9  dAb 


2nr 


(68) 


where  k  is  an  absorption  coefficient  and  r  is  the  distance  from  the  area  on  the  black 
body  surface  to  the  small  volume.  The  factor  in  the  parenthesizes  with  the  radius  of 
the  black  body  Rb  accounts  for  the  curvature  of  the  surface.  This  expression  is  to 
be  integrated  over  the  entire  planar  surface  of  the  black  body  by  using  cylindrical 
coordinates  (p,  <p).  The  area  dAb  is  then  p  dp  dcp. 


K{Tb’s) 

ne{Tb,s) 


NF(Tb ) 


2  n 


R, 


i  In  co 


\  Rb  +  5 , 


NF(Th) 


R, 


\Rb  +  s , 


W,  exp(-kr)  cos 9 

P  dp - 1 — f - 

0  0  r 

r  exp  (-kr)  cos9 

^  A - 


(69) 
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for  the  total  production  rate  for  the  photoelectrons  per  unit  volume.  The  integral  is 
evaluated  by  changing  the  variable  for  integration.  First  let  p  dp  =  r  dr  from  the 
differentiation  of  r2  =  p2  +  s2  and  cos  6  =  s/r, 


ex 


P  j~kr) 


(70) 


and  then  let  x  =  kr  for  dr  =  dx/k, 

I{ks)  =  ks]  dxexp(~x\ 

ks 


and 


ne(Tb,s)  =  N 


f  Rb  ' 

\Rb+s  J 


F{Th)l{ks). 


(71) 


(72) 


The  total  production  rate  for  the  region  between  the  black  body  surface  and  the 
distance  s  is 


he{s)=  N 


R, 


\Rb  +  s  j 


P(Tb)(\-l(ks)). 


(73) 


The  integral  I(ks)  was  numerically  integrated  as  a  function  of  ks,  which  is  the  optical 
density  of  the  layer.  The  values  of  the  integral  were  fitted  to  an  empirical  formula, 

I{ks)  =  0.6906  exp  (-1.5  2  As)  +  0.3094  exp  (-8.32As) ,  (74) 


that  is  a  sufficient  approximation  for  the  integral.  If  it  is  assumed  that  the  spectrum 
of  the  photons  is  proportional  to  the  black  body  throughout  the  thick  layer,  then  a 
weighted  cross  section  may  be  defined  by  first  substituting  Eq.  57  for  Np(s,  hv)  in 
Eq.  64  and  integrating  over  the  photon  energy, 


and  then  define 

oo  oo 

Qb(T»)\P(hv,Tt)  d(hv)=  Je(Av)  P(hv,Tt)  d(hv) 

hu0  huQ 

for  a  cross  section  that  is  weighted  over  the  photon  energy.  Also  let 


(75) 


(76) 
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(77) 


r(r6)=  \p(hv,Tb)d(hv) 

huQ 

for  the  total  number  of  photons  per  area  per  second  with  energies  greater  than  hv0. 
With  these  definitions,  the  equation  is  then 

ive,(rs)r(r4),  (78) 

OS 

where  a  coefficient  of  absorption  can  be  defined  as 

|7  Q(hv)  P(hv,Tb )  d(hv) 

k  =  NQb{Tb)  =  N -  (79) 

j,  P(hv,Tb)  d(hv) 

J  hua 

for  the  black  body  distribution.  The  weighed  cross  section  per  atom  is  shown  in 
Fig.  10  for  an  all  molecular  air  where  huo  =  12.07  eV  is  the  ionization  potential  of 
O2.  The  weighted  cross  section  for  all  atomic  air  has  huo  =  13.62  eV,  which  is  the 
ionization  potential  for  atomic  O. 


Fig.  10  Weighted  cross  section  per  atom;  molecular  air  ( - )  and  atomic  air  ( - ) 

Since  it  is  expected  that  most  of  the  photoionization  will  occur  when  the  air  is 
mostly  molecular,  the  weighted  cross  section  in  Megabarns  was  fitted  to  the 
following  polynomial  for  temperatures  greater  than  4000  K, 
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(80) 


Qh  {Tb)=  2.209  x  10  i87;4  -  3.398  x  10  I37;3 

+  1.647x10  87;2  - 1.32 1x  10  4r,+ 1.064 

and  Qb(Tb)  =  0.778  Mb  when  the  temperature  is  less  than  4000  K. 

In  the  finite  difference  method  for  the  electric  arc  there  are  a  number  of  layers. 
Since  these  layers  may  have  different  thicknesses  and  different  absorption 
coefficients  to  give  different  optical  densities,  determining  the  photon  flux  for  each 
layer  would  be  difficult  when  the  procedure  for  the  single  layer  is  applied  for  a 
number  of  layers.  Thus,  an  assumption  is  made  to  ease  the  calculation.  When  the 
photons  leave  the  absorbing  layer,  they  have  a  different  angular  distribution  when 
entering  the  next  layer.  Photons  having  an  angle  close  to  jt/2  are  lessened  by  cos  6, 
and  they  have  a  long  absorption  path.  Photons  with  zero  angle  are  not  lessened  by 
cos  9,  and  they  have  a  shorter  absorption  path.  Depending  on  the  optical  density  of 
the  layer,  this  changes  the  angular  distribution  of  the  photons  and  attenuates  the 
photons.  Layers  with  low  optical  densities  less  than  about  0.1,  for  example,  have 
angular  distributions  close  to  that  of  the  black  body  but  attenuated  by  about  10%. 
Layers  with  high  optical  densities  greater  than  about  1 .0,  for  example,  have  angular 
distributions  that  are  more  pronounced  at  small  angles  but  attenuated  by  about  60%. 
Also,  since  the  optical  density  depends  on  the  frequency  of  the  photons,  the 
spectrum  of  the  photons  after  passing  through  the  layer  may  no  longer  be  that  of  a 
black  body.  Still,  it  is  assumed  that  the  spectrum  of  these  photons  is  a  black  body 
but  attenuated,  and  the  angular  distribution  of  the  photons  is  the  same  as  the  black 
body  so  that  these  equations  can  be  applied  to  the  next  layer. 

As  an  example  of  the  change  in  the  electrical  conductivity  of  a  layer  due  to 
photoionization,  a  layer  just  outside  the  channel  was  chosen  from  the  calculations 
that  produced  the  results  shown  in  Figs.  2-5.  At  a  time  near  the  peak  current,  6.0  ps, 
a  layer  had  the  following  properties:  temperature  (T  =  2436  K),  density 
(p  =  1.31x1  O'3  g  cm'3),  and  width  (s  =  8. 17x  10"4  cm).  The  layer  next  to  it  was  inside 
the  channel  and  had  a  black  body  temperature  of  34,000  K.  The  time  step  for  the 
calculation  was  1 .68x  10'9  s.  The  total  number  density  of  atoms  and  ions  in  the  layer 
is  N=  (1.31x10~3)(6.022x1023)/14.43  =  5.43x1019  cm-3,  where  14.43  is  taken  to  be 
an  atomic  weight  of  air.  Before  this  layer  is  exposed  to  the  black  body  radiation, 
the  fractions  for  the  components  of  air  from  Eq.  32  gives  Ao  =  2.89x10"7, 
Ai  =  5. 93x1 0'15,  and  A2  =  1.41xl0"63.  From  these  fractions,  the  total  number  density 
for  the  neutral  atoms  even  when  they  are  in  molecules  is  no  =  5.43x1019  cm'3.  The 
number  density  for  the  electrons  is  ne  =  0.09  cm'3,  and  the  number  density  for  the 
singly  charged  ions  is  m  =  0.09  cm'3.  The  number  density  for  the  doubly  charged 
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ions  is  m  =  1.3><1  O'64  cm~3.  Thus,  with  the  mean  ionic  charge  Z  (Eq.  35)  of  1.0  and 
temperature  of  2436  K,  the  electrical  conductivity  of  the  layer  (Eq.  37),  is 
7.27x10"18  S/cm.  Now  let  this  layer  be  exposed  to  the  black  body  radiation  at  a 
temperature  of  34,000  K,  which  will  produce  photoelectrons  with  the  energy 
equivalent  to  this  temperature.  The  weighted  cross  section  for  photoionization 
(Eq.  80)  is  5.21xlCT18  cm2-  for  a  coefficient  of  absorption 
k  =  no  5.21  xlO-18  =  283  cm-1  and  an  optical  thickness  ks  =  k  4.49x10'4  =  0.231. 
The  total  absorption  in  the  layer  is  I(ks)  =  0.531,  and  the  production  rate  is 
F(34000)  =  1.784  x  108  s_1.  The  total  production  rate  for  the  photoelectrons  per 
volume  per  second  in  the  layer  is  then 

ne  (5)  =5.43x  1019  (1.784  x  108)(l  -  0.531)  =  4.54  x  1027  (81) 

when  the  radius  of  the  black  body  Rt  is  larger  than  the  thickness  of  the  layer.  After 
the  time  step  of  1 .68  x  1 0  9  s,  the  number  density  of  the  photo  electrons  is 
ne  =  4.54xl027  (1.68xl0~9)  for  ne  =  7.63x1018  cm-3.  The  number  density  of  the 
neutrals  is  no  =  5 .43  x  1 0l9—ne  for  no  =  4.67x1019  cm  3,  and  the  number  density  of 
singly  charged  ions  is  m  =  0.09  +  ne  for  m  =  7.63x1018  cm-3.  The  number  density 
for  the  doubly  charged  ions  does  not  change.  The  mean  ion  charge  is  still  Z  =  1 
because  the  number  density  for  the  doubly  charged  ions  is  very  small.  Using  these 
number  densities  and  an  electron  temperature  of  34,000  K,  the  conductivity  is 
117S/cm,  which  is  to  be  added  to  the  electrical  conductivity  before 
photoionization:  7.27x  10'18  S/cm.  In  other  parts  of  the  arc,  the  conductivities  before 
and  after  photoionization  are  comparable.  Thus,  the  total  of  the  conductivities  in 
this  case  is  117  S/cm,  while  a  typical  electrical  conductivity  in  the  channel  is 
200  S/cm. 

This  chosen  layer,  however,  is  also  being  photoionized  by  an  attenuated  flux  of 
photons  from  other  more-distant  layers  having  different  temperatures  where  the 
temperature  of  these  photoelectrons  are  set  to  be  equal  to  the  temperature  of  the 
distant  layer.  The  attenuation  of  the  photon  flux  from  a  distant  layer  is  determined 
by  following  the  attenuation  through  each  intervening  layer  having  an  optical 
thickness  and  a  value  for  I(ks).  This  attenuated  flux  is  used  for  an  additional 
production  rate  for  the  photoelectrons  in  the  layer.  The  total  production  rate  is 
multiplied  by  the  time  step  for  the  number  density  of  the  photoelectrons  that  are 
added  to  the  free  electrons  in  the  layer.  A  like  number  is  added  to  the  number 
density  of  the  singly  charged  ions  and  subtracted  from  the  number  density  of  atoms. 
The  number  density  for  the  doubly  charged  ions  is  not  changed.  Since  it  is  assumed 
that  the  spectrum  from  the  distant  layer  is  a  black  body,  the  temperature  of  these 
photoelectrons  is  equal  to  the  temperature  of  the  distant  layer.  An  average  electron 
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temperature  in  the  /th  layer  produced  by  all  the  distant  layers  (i^j)  may  be  defined 
by  a  weighted  average, 


T,  = 


(".I, 


(82) 


where  the  added  half  integer  for  the  indexes  of  these  intrinsic  properties  is  dropped. 
This  temperature  is  then  used  to  define  a  temperature  for  all  of  the  photoelectrons 
by  2  sums.  Let  the  first  summation  be 

S;=Sr'+7’(n,).A(,,  (83) 


which  is  a  sum  of  the  product  of  the  average  temperature  of  the  newly  produced 
photoelectrons  and  their  number.  The  second  summation  is  the  number  of  the  newly 
produced  photoelectrons, 

Dr  =  D"-'  +(«,),  A/,.  (84) 

Atn  is  the  time  step  and  (lVe).  is  the  total  production  rate  of  the  additional 
photoelectrons  by  all  the  other  elements.  Thus,  the  average  temperature  for  all 
photoelectrons  may  be  defined  as 

_  r*« 

7;.=-*-.  (85) 

'  d; 


During  this  time,  however,  all  of  the  photoelectrons  in  an  element  may  be 
recombining  with  the  ions.  When  an  electron  comes  close  to  an  ion,  the  electron 
could  lose  energy  when  a  photon  is  emitted  and  recombine  with  the  ion.  This 
sequence  is  the  time  reverse  of  photoionization.  The  cross  section  for 
recombination  for  hydrogen-like  atoms13  as  a  function  of  the  speed  of  the  electrons 

Ve  is 


a  = 


1.1x10 


(86) 


where  the  cross  section  is  in  square  centimeters.  Although  the  ions  are  not 
hydrogen,  the  upper  energy  levels  of  the  oxygen  and  nitrogen  ions,  where  most  of 
the  electrons  will  recombine,  are  hydrogen-like.  The  recombination  rate  per  unit 
volume  and  per  unit  time  for  electrons  with  a  velocity  ve  is 

Nc  =  nn{(jcve.  (87) 
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ne  is  the  number  density  of  the  electrons  and  n  /  is  the  number  density  of  the  singly 
charged  ions.  If  the  electrons  have  a  Maxwell  velocity  distribution  for  a  temperature 
f(ve,T),  the  weighted  recombination  rate  for  the  electrons  is 


Rc  =  l.lx\0  5 nen^dve 


oo 


f(ve,T)  _2.2xlO  uneiil 


(88) 


0 


V. 


e 


Thus,  at  the  end  of  the  time  step,  the  number  density  in  the  element  for  the  electrons, 
neutrals,  and  singly  charged  ions  are 


(89) 


'+(*.- K)*. 


where  the  index  for  the  element  is  dropped  and  the  recombination  with  the  doubly 
charged  ions  are  ignored.  The  number  densities  for  the  neutrals  and  the  singly 
charged  ions  must  be  greater  than  zero  but  less  than  the  total  number  density  N.  The 
number  density  for  the  electrons  can  be  greater  than  the  total  number  density  when 
there  are  doubly  charged  ions  present.  These  limits  are  enforced  at  the  end  of  each 
time  step. 

6.  Results 


The  finite  difference  equations  (Eqs.  8-16)  were  solved  using  Eq.  28  for  the  thermal 
conductivity  for  2  cases.  The  initial  condition  for  each  case  has  the  air  within  a 
radius  of  0.1  cm  with  a  temperature  of  15,000  K  at  ambient  density.  The  total 
conductivity  of  the  arc  (Fig.  11)  was  calculated  for  a  current  given  by  Eq.  21.  The 
first  case  has  no  production  of  the  photoelectrons  (no  PEs)  and  Eq.  50  for  the 
electrical  conductivity.  The  second  case  has  the  PEs  and  Eq.  50  for  the  electrical 
conductivity.  The  temperature  and  pressure  at  the  center  of  the  arc,  the  radius  of  the 
channel,  and  the  position  of  the  shock  wave  are  compared  with  Akram’s  results,7,8 
which  include  radiation  transfer.  Fig.  1 1  is  the  conductivity  of  the  arc  as  a  function 
of  time  when  there  is  photoionization  (PE)  and  when  there  is  no  photoionization 
(no  PE).  When  there  is  photoionization,  the  atoms  in  the  channel  quickly  become 
completely  ionized,  which  results  in  making  the  channel  more  conductive.  The  total 
conductivity  using  Eq.  36  and  no  PE  is  also  shown  to  simulate  Akram’s  results. 
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Fig.  11  Conductivity  of  the  arc:  PE  (—  ),  no  PE  ( - ),  and  Eq.  36  (  ) 

The  temperature  at  the  center  of  the  arc,  x  =  0  (Fig.  12)  tends  to  decrease  as  the 
electrical  conductivity  increases.  This  is  reasonable  because  an  increase  in  the 
conductivity  results  in  a  decrease  in  the  Joule  heating  for  a  given  current.  The 
temperature  by  Akram,7’8  also  in  Fig.  12,  follows  this  trend.  The  electrical 
conductivity  used  by  Akram  (Fig.  6)  is  generally  higher  than  Eq.  50,  and 
presumably  the  inclusion  of  radiation  transfer  by  Akram  would  also  increase  the 
conductivity.  In  comparison,  using  Eq.  50  with  no  production  of  photoelectrons 
resulted  in  a  high  core  temperature. 


Fig.  12  Core  temperature:  PE  ( - ),  no  PE  ( - ),  and  Akram7  ( - ) 

The  pressure  at  the  center  of  the  arc  (Fig.  13)  seems  to  be  insensitive  to  the 
photoelectrons.  In  the  case  where  there  are  no  photoelectrons,  the  pressure  agrees 
with  the  pressure  according  to  Akram. 
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Fig.  13  Core  pressure:  PE  ( - ),  no  PE  ( - ),  and  Akram7’8  (  ) 

Figure  14  shows  the  radius  of  the  channel  and  the  position  of  the  shock  wave  during 
the  current  pulse.  The  radius  of  the  channel  is  taken  as  the  radius  where  the 
temperature  is  half  of  the  temperature  at  the  very  center  (Fig.  2).  The  position  of 
the  shock  wave  is  where  there  is  the  maximum  pressure  (Fig.  3).  The  position  of 
the  shock  wave  and  the  radius  of  the  channel  are  again  close  to  Akram’ s  results 
when  there  are  no  photoelectrons.  These  radii  are  slightly  reduced  when 
photoelectrons  are  included  because  of  the  reduced  Joule  heating  of  the  channel  due 
to  its  high  conductivity. 


Time  (|xs) 

Fig.  14  Radius  of  the  channel  and  shock  wave:  PE  ( - ),  no  PE  ( - ),  and  Akram7,8  (  ) 
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Except  for  the  core  temperature,  when  there  is  no  photoionization,  the  present 
results  agree  with  Akram’s  results.  Using  the  lower  conductivity  and  omitting  the 
photoionization  that  gave  a  higher  core  temperature  for  these  seems  to  be 
comparable  to  the  higher  conductivity  and  the  radiation  transfer  in  Akram’s  results. 
When  the  photoionization  was  included,  the  core  temperature  was  reduced  and  the 
conductivity  of  the  arc  increased.  Photoionization,  however,  had  little  effect  on  the 
pressure  at  the  center,  the  radius  of  the  channel,  and  the  position  of  the  shock  wave. 


0.0  0.5  1.0  1.5  2.0  2.5  3.0 

Radius  (cm) 


Fig.  15  Ion  temperature  ( - ),  electron  temperature  (  ),  and  relative  density  ( - )  at  25  ps 

When  there  is  photoionization,  the  temperature  of  the  electrons  is  different  from  the 
temperature  of  the  ions  and  atoms.  Figure  15  is  a  profile  of  the  electron  temperature 
and  the  ion  temperature,  which  includes  the  atoms  as  well.  The  relative  density  of 
air  is  also  included  with  its  scale  on  the  right  side.  These  profiles  are  at  a  time  of 
25  |as  when  high-temperature  electrons  in  dense  air  are  evident.  The  number  density 
of  these  electrons,  however,  are  relatively  low,  but  the  conductivity  of  the  air  is 
sufficient  for  Joule  heating  to  be  noticeable  at  the  radius  of  about  1.7  cm. 

7.  Conclusion 


The  results  of  this  study  evaluated  quantities  that  are  accessible  to  measurement 
and  whether  they  can  be  used  to  infer  the  mechanisms  in  the  arc.  The  radius  of  the 
conducting  channel  can  be  measured  by  laser  interferogram.3  The  peak  pressure  of 
the  shock  wave  can  be  measured  and  its  position  can  be  sensed  by  pressure  sensors. 
These  quantities,  however,  are  nearly  identical  in  2  cases.  In  the  first  case,7  the 
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electrical  conductivity  was  generally  larger  than  the  one  given  by  Eq.  50  and 
included  radiation  transfer.  The  second  case  used  Eq.  50  for  the  electrical 
conductivity  and  did  not  include  radiation  transfer  and  no  photoionization  (no  PE). 
Thus,  these  measurements  may  not  distinguish  the  2  cases.  The  temperature  at  the 
very  center  is  different,  but  this  temperature  is  not  assessable  to  measurement.  The 
total  conductivity  of  the  arc,  however,  is  different  and  assessable  to  measurement 
that  can  distinguish  these  2  cases.  Since  the  total  conductivity  of  the  arc 
conductivity  is  even  higher  than  these  cases  when  photoionization  is  included  (PE), 
the  effects  of  photoionization  can  be  evaluated.  Detennining  the  high  conductivity 
of  an  arc  is  difficult  because  it  requires  careful  measurements  and  analysis  of  the 
current,  the  voltage  across  the  arc,  and  the  radius  of  the  arc. 

It  was  assumed  in  the  previous  studies  that  a  finite  element  radiates  as  a  black  body. 
This  assumption  was  used  here  so  that  the  present  results  can  be  compared  with 
previous  calculations.  A  number  of  assumptions  were  made  to  ease  the  calculation 
of  the  radiation  transfer  via  photoionization  while  other  processes  were  ignored. 
This  assumption  resulted  in  different  temperatures  for  the  ions  and  electrons. 
Eventually,  these  temperatures  should  tend  to  become  equal  to  reach  thermal 
equilibrium.  Thus  a  more  detailed  description  of  the  radiation  transfer  and  the  other 
processes  will  be  included  in  the  future. 


Approved  for  public  release;  distribution  is  unlimited. 

32 


8.  References 


1 .  Engel  TG,  Donaldson  AL,  Kristiansen  M.  The  pulsed  discharge  arc  resistance 
and  its  functional  behavior.  IEEE  Transactions  on  Plasma  Science. 
1989;17(2):323-329. 

2.  Akiyama  H,  Kristiansen  M,  Krompholz  H,  Maas  B.  Current-voltage 
characteristics  of  a  high-current  pulsed  discharge  in  air.  IEEE  Transactions  on 
Plasma  Science.  1988;  16(2):3 12-3 16. 

3.  Kushner  MJ,  Kimura  WD,  Byron  SR.  Arc  resistance  of  laser-triggered  spark 
gaps.  J  Appl  Phys.  1985;58(5):  1 744—175 1 . 

4.  Plooster  MN.  Numerical  simulation  of  spark  discharges  in  air.  The  Physics  of 
Fluids.  1971;14(16):21 11-2123. 

5.  Plooster  MN.  Shock  waves  from  line  sources:  numerical  solutions  and 
experimental  measurements.  The  Physics  of  Fluids.  1970;  13(1 1):2665-2675. 

6.  VonNeumann  J,  Richtmyer  RD.  A  method  for  the  numerical  calculation  of 
hydrodynamic  shocks.  J  Appl  Phys.  1950;21:232-237. 

7.  Akram  M,  Lundgren  E.  The  evolution  of  spark  discharges  in  gases:  I. 
macroscopic  models.  J  Phys  D  Appl  Phys.  1996;29:2129-2136. 

8.  Akram  M.  The  evolution  of  spark  discharges  in  gases:  II.  numerical  solution 
of  one-dimensional  models.  J  Phys  D  Appl  Phys.  1996;29:2137-2146. 

9.  Olsen  HN.  Thermal  and  electrical  properties  of  an  argon  plasma.  Physics  of 
Fluids.  1959;2(6):614-623. 

10.  Spitzer  L,  Harm  R.  Transport  phenomena  in  a  completely  ionized  gas.  Physical 
Review.  1953;89(5):977-981. 

1 1 .  Riviere  PH,  Soufiani  A,  Perrin  MY,  Riad  H,  Gleizes  A.  Air  mixture  radiative 
property  modelling  in  the  temperature  range  10,000-40,000  K.  J  Quant 
Spectrosc  Radiat  Transfer.  1996;56(l):29-45. 

12.  Chauverau  S,  Deron  C,  Perrin  M,  Riviere  P,  Soufiani  A.  Radiative  transfer  in 
LTE  air  plasmas  for  temperatures  up  to  15,000  K.  J  Quant  Spectrosc  Radiat 
Transfer.  2003;77:113-130. 

13.  Zeldovich  YB,  Raiser  YP.  Physics  of  shock  waves  and  high-temperature 
hydrodynamic  phenomena.  In:  Hayes  WD,  Probstein  RF,  editors.  Temperature 
hydrodynamic  phenomena.  Mineola  (NY):  Dover  Publications,  Inc.;  2002. 


Approved  for  public  release;  distribution  is  unlimited. 


33 


14.  Press  WH,  Flannery  BP,  Teukolsky  SA,  Vetterling  WT.  Numerical  recipes. 
(London):  Cambridge  University  Press;  1986. 

15.  Plooster  MN.  Shock  waves  from  line  sources.  Boulder  (CO):  National  Center 
for  Atmospheric  Research;  1968  Nov.  Report  No.:  NCAR-TN-37. 

16.  Richtmyer  RD,  Morton  KW.  Difference  methods  for  initial  value  problems. 
2nd  ed.  New  York  (NY):  Interscience;  1994.  p.  405. 

17.  Lamb  L,  Lin  S.  Electrical  conductivity  of  thermally  ionized  air  produced  in  a 
shock  tube.  J  Appl  Phys.  1957;28(7):754. 

18.  Hirshman  SP.  Transport  of  a  multiple  ion  species  plasma  in  the  Pfirsch- 
Schluter  regime.  Physics  of  Fluids.  1977;20(4):589-598. 

19.  Callen  JD.  Fundamentals  of  plasma  physics.  Draft;  2003  [accessed  2016  July 
26] .  http://homepages.cae.wisc.edu/~callen/book.html. 

20.  Cohen  RS,  Spitzer  L,  Routly  PM.  The  electrical  conductivity  of  an  ionized  gas. 
Physical  Review.  1950;80(2):230-238. 

21.  Itikawa  Y,  Hayashi  A,  Ichimura  K,  Onda  K,  Sakimoto  K,  Takayanagi  K, 
Nakamura  M,  Nishimura  H,  Takayanagi  T.  Cross  sections  for  electrons  and 
photons  with  nitrogen  molecules.  J  Phys  Chem  Ref  Data.  1986;  15(3):985— 
1010. 

22.  Itikawa  Y,  Ichimura  A,  Onda  K,  Sakimoto  K,  Takayanagi  K,  Hatano  Y, 
Hayashi  M,  Nishimura  H,  Tsurubuchi  S.  Cross  sections  for  collisions  of 
electrons  and  photons  with  oxygen  molecules.  J  Phys  Chem  Ref  Data. 
1 989;  1 8(  1):23— 42. 

23.  Samson  JAR,  Angel  GC.  Single-  and  double-photoionization  cross  sections  of 
atomic  nitrogen  from  threshold  to  3 1  A.  Physical  Review  A.  1990;42(1):  1307— 
1312. 

24.  Angel  GC,  Samson  JAR.  Total  photoionization  cross  sections  of  atomic 
oxygen  from  threshold  to  44.3  A.  Physical  Review  A.  1988;38(1 1):5578- 
5585. 


Approved  for  public  release;  distribution  is  unlimited. 

34 


List  of  Symbols,  Abbreviations,  and  Acronyms 


1-D 

1 -dimensional 

eV 

electron  volt 

Mb 

Megabarn 

N 

nitrogen  atom 

n2 

nitrogen  molecule 

no  PE 

no  photoelectrons  produced 

0 

oxygen  atom 

02 

oxygen  molecule 

PE 

photoelectrons  produced 
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